Growth Rate of Large Scale Structure as a Powerful Probe of Dark Energy 



cn ■ 
O ■ 
o : 

Oh 
< 



> ■ 

oo ; 

^ . 

^ ■ 

o : 
m 
o • 

Oh; 

6 



Asantha Cooray/ Dragan Huterer,^ and Daniel Baumann'^ 

^Theoretical Astrophysics, California Institute of Technology, Pasadena, CA 91125 
^Department of Physics, Case Western Reserve University, Cleveland, OH 44^06 
^ DAMTP, Centre for Mathematical Sciences, University of Cambridge, Cambridge CB3 OWA, UK 

The redshift evolution of the growth rate of the gravitational potential, d[D / a) / dz, is an excellent 
discriminator of dark energy parameters and, in principle, more powerful than standard classical 
tests of cosmology. This evolution is directly observable through the integrated Sachs- Wolfe effect 
in cosmic microwave background (CMB) anisotropies. We consider the prospects of measuring the 
growth rate via a novel method employed through measurements of CMB polarization towards 
galaxy clusters. The potentially achievable errors on dark energy parameters are comparable and 
fully complementary to those expected from other upcoming tests of dark energy, making this test 
a highly promising tool of precision cosmology. 



One of the key issues in modern cosmology is devel- 
oping efficient and complementary methods to measure 
cosmological parameters and cosmological functions. In 
particular, much interest has been devoted to developing 
methods to constrain the properties of the mysterious 
dark energy component that causes the recently discov- 
ered accelerated expansion of the universe P]. To this 
end, it has been pointed out that type la supernovae, 
number counts, and weak gravitational lensing are all 
very promising probes of the dark energy equation of 
state w and its energy density relative to critical J^de 
(see and references therein), and that a number of 
other methods are likely to contribute useful information. 

These cosmological tests probe various fundamental 
quantities. For example, type la supernovae effectively 
measure the luminosity distance, number counts are sen- 
sitive to a combination of the volume element and the 
growth of density perturbations, while cosmic microwave 
background (CMB) anisotropy effectively determines the 
distance to the surface of last scattering. While these 
tests are well understood and pursued in various obser- 
vational programs, the information one can extract from 
these measurements is limited by the presence of fun- 
damental degeneracies of cosmological parameters that 
enter the observable quantity in question. 

It is typically advantageous when the measurements 
involve not the quantity itself but rather its derivatives 
with respect to time or redshift, since in that case the 
dependence on the equation of state ^(z) is more direct. 
For example, the Hubble parameter H{z) is more sensi- 
tive to the equation of state than the comoving distance 
r{z) since r{z) — J dz/H{z), but the latter has the ad- 
vantage of being readily and accurately measurable. In 
this respect, the linear growth factor of density pertur- 
bations provides important information since it is a func- 
tion of the Hubble parameter and the equation of state 
of dark energy w (see below). 

As we will show, the rate of evolution of the growth 
factor with redshift is a tremendous tool for measur- 
ing dark energy parameters. We will further suggest 
the integrated Sachs- Wolfe (ISW; 0) effect as a probe 
for this purpose and propose polarization measurements 



of CMB anisotropy towards galaxy clusters. The latter 
provides an indirect method to extract the temperature 
quadrupole associated with the ISW effect as a function 
of the cluster redshift, with a reduction in the cosmic vari- 
ance which plagues large scale temperature anisotropy 
measurements. To make this study practical, we con- 
sider the prospects of upcoming arcminute scale CMB 
polarization observations with instruments such as the 
South Pole Telescope (SPT) and the planned CMBPol 
satellite mission. 

To begin, we review aspects related to the growth of 
large scale structure. In linear theory, all Fourier modes 
of the density perturbation, d{= Spu/PrmM), grow at 
the same rate: Sk{a) = D{a)Sk{a — 1), where D{a) is 
the growth factor normalized to unity today and a = 
(l-|-z)~^ is the scale factor. In the matter-dominated era 
D{a) — a, while in the presence of a smooth dark energy 
component perturbation growth slows and D{a) increases 
less rapidly with a. In general, the growth function can 
be computed by solving the linear perturbation equation 
Sk + 2{a/a)5k — 'inGpuSk = where dot is the derivative 
with respect to physical time. 

Defining the growth suppression rate (growth rate rel- 
ative to that in a flat, matter-dominated universe) as 
g{a) = D{a)/a, and still allowing for a general w(a), one 
can write 



din 



[5 - 3w{a)flDE{a)] - 
a ma 

3 [1 — w{a)] fluEia)g — 0, 



(1) 



where r2DE(a) is the fractional dark energy density at the 
scale factor a. For constant w, the solution is given in 
terms of the hypergeometric function Q, while to com- 
pute g{a) and/or D{a) for a non-constant w{a) one can 
either solve Eq. numerically or use analytic approxi- 
mations 

It has long been known that the growth function 
strongly depends on f^M, the fractional density in matter, 
and w. Also, the strength of several cosmological tests, 
such as number counts |y| , clustering measured in redshift 
slices and weak lensing comes primarily from their 
dependence on the growth function D{z). On the other 
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FIG. 1: Errors in the JIm-io plane (top panel) and the uiQ-iOa 
plane (bottom panel, with a prior on matter density, G(Q.y\) 
of 0.01) assuming a 10% measurement of a given quantity at 
redshifts z = [0.1, 0.2, . . . , 2.0]. We show errors for distance 
r(z), Hubble parameter H{z), growth factor D{z), differential 
volume element r^(z)/H{z), dD/drj [t] is the conformal time), 
and rate of growth suppression, dg/dz. Note that the latter, 
dg/dz, is by far the most sensitive to and w. For exam- 
ple, a 10% accuracy measurement of dg/dz is comparable and 
complementary to < 1% measurement of distance. 



hand, it has been known that redshift or time deriva- 
tives of distance are more directly related to dark energy 
parameters; in particular, the equation of state w{z) is 
directly related to the first and second derivatives of dis- 
tance with respect to redshift j^]. Unfortunately, the 
derivatives are not directly measured but are obtained 
by taking numerical derivatives of noisy data, which sig- 
nificantly increases the error in the reconstructed w{z). 

It is interesting to examine the sensitivity of the rate 
of change of the growth suppression factor to f^M and 
w. To do this, we first consider constant w, and then a 
two-parameter description of time-varying w. For the 
latter we do not choose the commonly used w{z) = 
Wo + w'z 0, but rather w{z) = wq + Waz/{\ + z) ^ 
which is bounded at high redshift and facilitates the inte- 
gration of Eq. (Q. [For w' aficionados, we mention that 



the error in Wa is roughly twice the error in w' ] Fig. ^ 
shows the error bars in the ^yi~w plane (top) and wq- 
Wa plane (bottom) using various classical tests assum- 
ing a fiducial model of flu = 0.3 and w = —1; we use 
the same fiducial model throughout the paper. The cal- 
culation uses the Fisher matrix formalism and assumes 
10% measurements in a given quantity at each interval of 
Az = 0.1 in redshift between z = 0.1 and z = 2. We show 
a variety of quantities, including the distance r(z), the 
volume element r'^{z)/H{z) and the growth factor D{z), 
which are the most commonly considered probes of dark 
energy. We also consider dD/drj ( rj i s conformal time) 
and dg/dz. As emphasized in Ref. [iJI, dD/drj, which is 
measured by large-scale velocities, is mostly sensitive to 
ilyi and not w. What Fig. ^ illustrates is that dg/dz is 
much more powerful than other probes due to the spe- 
cific way the degeneracy is broken. For example, for the 
same relative accuracy in observations, dg/dz is about 
15 times stronger than the comoving distance r{z)\ Of 
course, this comparison is not necessarily fair, since enor- 
mous amount of work has gone into developing methods 
to determine distances, which are now expected to be 
measurable to an accuracy of about 1% (per interval of 
0.1 in redshift) by SNAP [IJj, making them the most di- 
rect probes of the cosmological expansion history, while 
not much attention has been devoted to the more esoteric 
quantity dg/dz. In the remainder of this paper we show 
that there indeed exists a very promising cosmological 
test which is sensitive to dg/dz. 

The above discussion indicates that it would be ideal to 
have a cosmological probe of the evolution of growth sup- 
pression, dg/dz. It turns out that just such a probe exists 
in a universe that is not matter-dominated at late times. 
The dark energy domination causes the time-variation 
of the gravitational potential, which in turn contributes 
to CMB anisotropics through the ISW effect 0. The 
resulting temperature fluctuation is given by 
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where Tree is the radial comoving distance to last scat- 
tering with Zrec = 1100. From Poisson's equation, 
V^<I> = 3/2 HgUuiS/a), it follows that the gravitational 
potential $ is proportional to the growth suppression g. 
The ISW effect therefore gives a direct measure of the in- 
tegral of dg/dr (or dg/dz) computed over some effective 
time (or redshift) interval. 

While the ISW effect determined at t he p resent time 
can be used as a probe of dark energy 14], its contri- 
bution to CMB temperature fluctuations is dwarfed by 
the primordial anisotropy contribution at last scatter- 
ing. Though the cross-correlation between the large scale 
structure and CMB anisotropy fluctuations has been con- 
sidered as a method to extract the ISW contribution |0| , 
such correlations are affected by the dominant noise con- 
tribution related to primary anisotropics [T^ . 

There is another way of extracting information cap- 
tured in the ISW effect: through the measurement of 
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CMB polarization towards galaxy clusters ^3 ■ The po- 
larization signal is generated by rscattering of the tem- 
perature quadrupole seen by free electrons in the cluster 
frame ^T§\. Provided that the optical depth to scatter- 
ing in individual clusters is determined a priori by other 
methods, such as the Sunyaev-Zel'dovich (SZ; [1^) ef- 
fect, one can measure the quadrupole at the cluster loca- 
tion with a reduction in cosmic variance ^^20] ■ Note that 
the quadrupole measured from a cluster at high redshift 
is not the same quadrupole as one observes today due 
to the difference in the projected length scales. Since 
the ISW effect contributes a significant fraction of the 
quadrupolar anisotropy at late times, cluster polariza- 
tion provides an indirect probe of dark energy. Because 
clusters can be selected over a wide range in redshift, 
the polarization signal can be measured as a function of 
redshift and inverted to reconstruct the evolution of the 
ISW quadrupole . 

The anisotropy quadrupole, C2{z), has two contribu- 
tions: one at the surface of last scattering due to the 
Sachs- Wolfe (SW) effect, CfW(z), and another at late 
times due to the ISW effect, Cl^^{z). We write these 
two contributions to the power spectrum, projected to a 
redshift z, respectively as 

CT'i^) ^ yI T^**(^'''-)j|[fc(rrcc-r)] and 
f°° dk 

Cf"(z) = Kiir / -^A|,(i:,r„e)x (3) 

Here r is the radial comoving distance out to redshift z 
and A|^(fc, ri.ec)(= A:'^P$$(fc, rrec)/27r^) is the logarith- 
mic power spectrum of fluctuations in the potential field 
at the last scattering surface. We will concentrate on 
the dark energy properties, whose effects are dominant 
at low redshifts, and assume that the parameters that de- 
fine the power spectrum, such as the normalization, spec- 
tral tilt, and physical matter and baryon densities flmh'^ 
and ile^^, are known to the accuracy expected from the 
Planck mission with polarization information .21.] . Given 
these priors, the SW contribution is then known to a few 
percent accuracy. Also note that, conveniently, only the 
large scales in the power spectrum contribute to the ISW 
effect, so that we do not need to consider thorny issues re- 
lated to small-scale non-linear structures and additional 
parameters such as the neutrino mass. 

The galaxy cluster polarization signal arises from the 
rescattering of the quadrupole which receives a contribu- 
tion from Cf^{z) at low redshifts. Ref. [13 discussed 
how well this quadrupole can be measured as a function 
of redshift with Planck and a ground-based experiment 
with significant instrumental noise. In the top panel of 
Fig. [5] we show the projected ISW contribution to the 
temperature quadrupole as a function of redshift, and 
expected errors for a ground-based survey targeting clus- 
ters down to a mass limit of IO^^Mq in a total area 
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FIG. 2: Top panel: the projected ISW quadrupole as a func- 
tion of redshift. The solid error bars assume a reconstruc- 
tion with clusters down to 10^^ Mq in an area of 10^ deg^ 
with an instrumental noise of 0.1 fiK. The dotted lines show 
the cosmic-variance for an all-sky reconstruction computed 
from the number of independent volumes sampled by clus- 
ters at each redshift bin |2C| . Bottom left: Parameter errors 
from the projected ISW quadrupole measurements, assuming 
w — const. The small ellipse is for the case shown in top 
panel (with cosmic variance added in quadrature), while the 
two larger ellipses assume a factor of 3 and 10 increase in 
the instrumental noise contribution, respectively. For com- 
parison, we also show the constraints expected from SNAP. 
When the most optimistic polarization information is added, 
snap's constraints on w improve by a factor of 3. Bottom 
right: same, but for wo and Wa and assuming an additional 
prior a{^lM) ~ 0.01. 



of 10^ deg^ with an instrumental noise for polarization 
observations of 0.1 /iK. As in 117|, we assume four chan- 
nels for these observations so that the ISW quadrupole 
can be separated from the contribution of the kinematic 
quadrupole. The latter has a distinct spectrum and the 
separation based on frequency information leads to an 
overall increase in noise by a factor of 2 to 3 depending on 
the exact frequencies of channels selected. Note that we 
have assumed an instrumental noise of 0.1 iiK for these 
observations. While a polarization sensitive detector ar- 
ray on the SPT can be expected to reach noise levels of ~ 
1 fiK or less per pixel, we have assumed an order of mag- 
nitude reduction in noise, as expected from the planned 
CMBPol satellite mission. Since the expected noise level 
for arcminute scale polarization observations from such a 
mission is not currently defined, and to consider ground- 
based efforts such as the SPT, we have considered the 
range of values between 0.1 and 1 so as to obtain 
some guidance on how well cluster polarization measure- 
ments with noise in this range can be used to probe dark 
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energy. 

In addition to instrumental noise, the polarization 
measurements are subject to cosmic variance. This vari- 
ance is determined by the number of independent vol- 
umes that last scattering spheres of individual clusters, 
in some redshift bin, occupy [23|- Dotted lines in the 
top panel of Fig. |21 show the cosmic variance contribu- 
tion for an all-sky experiment. As one moves to higher 
redshift, the number of independent samplings of the lo- 
cal quadrupole increases, leading to a reduction in cosmic 
variance. The expected redshift distribution of clusters 
peaks at redshifts around 1-1.5 where it provides the best 
estimate of the local quadrupole, while errors increase at 
very low and high redshift due to the smaller number of 
clusters. 

To consider how well these observations can be used to 
understand dark energy parameters, we again perform a 
Fisher matrix calculation. The bottom panels of Fig. |21 
show how well J^m and w (assuming a flat universe and 
constant w), and wq and Wa (assuming a two-parameter 
description of w{z) as before and a prior on flyi of 0.01) 
can be measured. While the errors are fairly large with a 
1 /zK noise level per pixel, improving this noise threshold 
to 0.1 iiK leads to significant gains in the determination 
of and w. Note also that these errors roughly scale 
as the inverse square root of the area of sky covered, 
and with all-sky coverage the errors are expected to de- 
crease by a factor of two. With an order of magnitude 
improvement in noise, the redshift evolution of the ISW 
effect extracted from polarization measurements becomes 



a powerful probe of dark energy providing significant es- 
timates of parameters, comparable and complementary 
to type la supernovae. 

To conclude, we have argued that the rate of evolu- 
tion of the growth suppression factor, dg/dz, is a very 
powerful probe of dark energy. We have shown that the 
polarization signal from a large number of galaxy clus- 
ters is directly related to this quantity, and can be used 
to constrain dark energy parameters. In the next decade, 
the planned mission CMBPol is expected to reach a sen- 
sitivity of order 0.1 /xK at arcminute resolution and have 
all-sky coverage, providing polarization measurements of 
a significant number (~ 10^) of clusters, from which the 
quadrupole can be reconstructed as a function of red- 
shift. Although our study is preliminary, we have shown 
that this method can provide constraints on the dark 
energy equation of state and its time variation compara- 
ble and complementary to those from type la supernovae 
and other well-studied probes of dark energy. More im- 
portantly, this method is entirely different from most of 
the others both in its theoretical underpinnings and in 
the systematic errors expected. Combining this method 
with others opens the exciting possibility of significantly 
improving the constraints on w and helps usher a new 
era in our exploration of dark energy. 
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